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The coalescence of a stellar-mass compact object into an intermediate-mass black hole (inter- 
mediate mass-ratio coalescence; IMRAC) is an important astrophysical source for ground-based 
gravitational- wave interferometers in the so-called advanced (or second- generation) configuration. 
However, the ability to carry out effective matched-filter based searches for these systems is limited 
by the lack of reliable waveforms. Here we consider binaries in which the intermediate-mass black 
hole has mass in the range 24 — 200 with a stellar-mass companion having masses in the 
range 1.4 — 18.5 M©. In addition, we constrain the mass ratios, q, of the binaries to be in the 
range 1/140 < ^ < 1/10 and we restrict our study to the case of circular binaries with non-spinning 
components. We investigate the relative contribution to the signal-to-noise ratio (SNR) of the three 
different phases of the coalescence: inspiral, merger and ringdown. We show that merger and ring- 
down contribute to a substantial fraction of the total SNR over a large portion of the mass parameter 
space, although in a limited portion the SNR is dominated by the inspiral phase. We further iden- 
tify three regions in the IMRAC mass-space in which: (i) inspiral-only searches could be performed 
with losses in detection rates L in the range 10% ^ L ^ 27%, (ii) searches based on inspiral-only 
templates lead to a loss in detection rates in the range 27% < L < 50%, and (iii) templates that 
include merger and ringdown are essential to prevent losses in detection rates greater than 50%. 
We investigate the effectiveness with which the inspiral-only portion of the IMRAC waveform space 
is covered by comparing several existing waveform families in this regime. We find that different 
waveform families are only marginally effective at describing one another, as measured by the "fit- 
ting factor". Our results reinforce the importance of extensive numerical relativity simulations of 
IMRACs to validate and calibrate semi-analytical waveform families and the need for further studies 
of suitable approximation schemes in this mass range. 

PACS numbers: 



I. INTRODUCTION 

Observations of ultra-luminous X-ray sources and sim- 
ulations of globular cluster dynamics suggest the exis- 
tence of intermediate-mass black holes (IMBHs). How- 
ever, observational evidence for their existence is still un- 
der debate, see e.g. [1, 2 . Gravitational waves from bi- 
nary coalescences involving IMBHs with masses 50 Mq 
^ M ^ 500 Mq are potentially detectable by advanced 
detectors - including Advanced LIGO [3j, Advanced 
Virgo U , and K AGRA [5] - with a low frequency cutoff 
of around 10 Hz. If IMBHs do exist, one likely contri- 
bution to gravitational-wave detections is believed to be 
through the coalescence of a compact stellar-mass com- 
panion (black hole or neutron star) with an IMBH, at 
a possible rate of up to -lOyr-i [Ms]. We will de- 
note these signals as intermediate mass-ratio coalescences 
(IMRACs) 1. 

Given that IMBHs in this mass range have proved ex- 



In the literature, the term frequently used for this class of ob- 
jects is intermediate mass-ratio inspirals or IMRIs, see e.g. [71 18]. 
However, in the context of ground-based observations, in partic- 
ular with second-generation instruments, we will show that the 
full coalescence is important for these systems, and it therefore 
seems to be more appropriate to call them IMRAC. 



tremely difficult to observe in the electromagnetic spec- 
trum, gravitational-wave detections may provide the first 
unambiguous observations of such objects through the 
robust measurement of their masses. Such observations 
would form an important channel for probing the dy- 
namical history of globular clusters. Furthermore, Ad- 
vanced LIGO/Virgo (aLIGO/ Virgo) may be able to pro- 
vide measurements of the quadrupole moment of a black 
hole [7, 9 , which would allow a null- hypothesis test of 
the Kerr metric for IMBHs. 

The gravitational waveform from the coalescence of 
two compact objects can be divided into three phases: 
a gradual inspiral, a rapid merger, and the quasi-normal 
ringdown of the resulting black hole. The relative contri- 
bution to the expected coalescence signal from inspiral, 
merger and ringdown is an important consideration for 
gravitational- wave searches. To leading Newtonian order 
the gravitational wave frequency at the inner-most stable 
circular orbit (ISCO) is 

/isco=4.4kHz (^^) , (1.1) 

where M is the total mass of the binary. For advanced 
detectors with a low frequency cut-off of ~ 10 Hz, we 
may only have access to either the very late stages of the 
inspiral, or solely merger and ringdown for the heaviest 
IMRAC systems. While the power in the merger and 
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ringdown is suppressed by a factor of the mass ratio rel- 
ative to the power in the inspiral, the fact that IMRAC 
systems are hable to merge either in-band, or at the low 
frequency limit of the bandwidth, means that merger and 
ringdown may be significant over a large portion of the 
detectable mass- space. Additionally, for cases where IM- 
RAC waveforms are inspiral-dominated, it is useful to 
know where inspiral-only searches could be targeted. 

Detecting IMRACs through gravitational waves 
will require template gravitational- wave form families 
adapted to highly asymmetrical mass-ratio systems. 
However the development of numerical relativity simu- 
lations and perturbative techniques in this regime is at 
an early stage which is potentially problematic (see, e.g. 
[? ] for a discussion of this issue). The issue of appro- 
priate template waveform families is thus central to the 
detection of IMRACs through gravitational waves. 

The effective-one-body approach, calibrated to numer- 
ical relativity, has led to template waveforms, known as 
EOBNR [10 , that describe the full inspiral, merger and 
ringdown coalescence-signal for comparable mass-ratio 
binaries; EOBNR waveforms should also be accurate at 
extreme mass ratios. However, to date only one full nu- 
merical simulation exists for mass-ratio q = 1/100 bi- 
naries [11]. Furthermore, EOBNR waveforms have been 
constructed to reproduce the dynamical evolution of bi- 
naries with extreme mass-ratios. EOBNR waveforms 
have not yet been compared to numerical relativity sim- 
ulations at such mass ratios, so their validity in the IM- 
RAC regime remains to be demonstrated. 

Meanwhile, in the context of extreme mass-ratio bina- 
ries, several authors have modelled the two-body motion 
by computing radiative and conservative self-force cor- 
rections to Kerr geodesic motion [121 US]- Waveforms 
computed within this scheme are inspiral-only and are 
only developed to lowest order in the mass ratio. These 
waveforms have been adapted to describe intermediate 
mass-ratio inspirals by including higher-order-in-mass- 
ratio corrections in [14 and have been used to study 
the detection of intermediate mass-ratio inspirals in the 
context of the proposed third-generation ground-based 
gravitational-wave interferometer the Einstein Telescope 
[15 . We refer to these intermediate mass-ratio inspiral 
waveforms as the "Huerta Gair" (HG) waveform family 
after its authors. This waveform family should be physi- 
cally well motivated to describe the inspiral of IMRACs. 

Typically one does not have an exact representation of 
"true" gravitational-wave signals but requires templates 
which are sufficiently effective at filtering such signals. 
A common metric for quantifying how well approximate 
waveform families are at filtering gravitational-wave sig- 
nals is known the "effectiveness", or fitting factor [16\ . 
This measures the fraction of the theoretical maximum 
signal-to- noise ratio (SNR) that could be recovered by 
using non-exact template waveforms. 

The work in this paper proceeds as follows. Firstly, 
by computing the effectiveness of inspiral-only template 
waveforms at filtering the full coalescence signal, we 



determine the relative importance of the inspiral and 
merger-ringdown phases. We identify three regions in 
the component mass plane in which: (a) inspiral-only 
searches are feasible with losses in detection rates L in 
the range 10% ^ L ^ 27%, (b) searches are limited by the 
lack of merger and ringdown in template waveforms and 
are liable to incur losses in detection rates in the range 
27% ^ L ^ 50%, and (c) merger and ringdown are es- 
sential for searches in order to prevent losses in detection 
rates greater than 50%. 

Secondly, to gain insight into the accuracy of the in- 
spiral portion of IMRAC waveforms we compute the ef- 
fectiveness of the inspiral-only portion of EOBNR wave- 
forms at filtering gravitational- wave signals as described 
by the HG waveform family. We find that there is a 
non-negligible discrepancy between EOBNR and HG in- 
spirals in the regime where inspiral-only searches could 
be considered sufficient. For reference we also compare 
EOBNR inspirals to a post-Newtonian (PN) waveform 
family known as TaylorT4 [16|. The PN expansion is li- 
able to be a poor choice of approximant for IMRACs be- 
cause of the large number of cycles spent at small radii. 
We find that EOBNR and HG are in better agreement 
with each other than to TaylorT4, as might be expected 
from the previous observation. 

Our approach does not directly address the accuracy of 
template waveforms, because none of the waveforms con- 
sidered have been matched to full numerical waveforms. 
However, assuming that the waveform families we con- 
sider "bracket" the correct gravitational waveforms in the 
IMRAC regime, this approach provides a useful heuristic 
for quantifying the effectiveness of existing gravitational 
waveforms for IMRAC searches. Further numerical rela- 
tivity simulations will be important in the continuing de- 
velopment of accurate template waveforms for IMRACs. 

Our analysis improves upon previous work to deter- 
mine the detectability of IMRAC sources [17] which only 
considered the so-called "faithfulness" of template wave- 
forms, i.e., the effectiveness of template waveforms eval- 
uated at the signal parameters. Additionally, that study 
only considered inspiral-only waveforms and focused on 
low frequency observations, e.g. with the proposed Laser 
Interferometer Space Antenna (LISA). 

This paper is organized as follows. In Sec. |TT] we de- 
scribe the waveform families used in our study. In Sec. [HI 
we compute the contributions to the SNR from the in- 
spiral and merger and ringdown phases of EOBNR wave- 
forms in the intermediate mass-ratio regime. In Sec. IV 
we study the effectiveness of inspiral-only waveforms to 
filter full coalescence signals from IMRAC sources and 
identify the three regions in which different searches 
could be conducted. In Sec. |V| we compare the inspi- 
ral portion of EOBNR waveforms to HG and TaylorT4 



waveforms. In Sec VI we consider the implications of our 
results for future searches in advanced detectors. 



3 



II. WAVEFORMS 

In this section we summarise the key concepts enter- 
ing the construction of the waveforms used in this study. 
Throughout the paper, for a binary system with indi- 
vidual component masses mi and m2 (with m2 < mi) 
we define the total mass as M = mi + m2, and mass 
ratio and symmetric mass ratio as = 7x12/ mi and 
rj = 17111712/ {mi + m2)^, respectively. 

We consider the family of waveforms constructed by 
calibrating the effective-one-body approach to numerical 
relativity (EOBNR) [10 . The EOBNR family describes 
the full inspiral-merger-ringdown signal; it is currently 
used in searches that reach the IMBH mass range, so 
far up to 100 Mq [18]. The free parameters in the fam- 
ily have been fitted to comparable mass ratio numerical 
relativity simulations, and by construction this family is 
deemed to be faithful in the test particle limit. For this 
work, we use the implementation provided by the LIGO 
Scientific Collaboration Algorithm Library (LAL) that 
corresponds to the approximant E0BNRv2. 

We also consider a waveform family based on test parti- 
cle motion in Kerr/Schwarzschild space-time with radia- 
tive and conservative self-force corrections which we refer 
to as the Huerta-Gair (HG) family [15 . The approxima- 
tion scheme is constructed specifically to handle highly- 
asymmetrical mass-ratio binaries and is therefore a phys- 
ically well motivated approximation scheme for interme- 
diate mass-ratio inspirals. These waveforms have been 
compared against, Teukolsky based waveforms for inspi- 
ralling test particles on geodesic orbits and the match 
exceeds 95% over a large portion of the parameter space 
[14]. These waveforms have been used to study detec- 
tion of intermediate mass-ratio inspirals in Einstein tele- 
scope in [15 . The Huerta-Gair waveforms describe only 
the inspiral portion of the coalescence signal. There is 
no corresponding LAL approximant. The gravitational- 
wave polarization states can be computed from Eqs. (14) 
and (15) of [15]. For our study, effects of orientation of 
the gravitational-wave source are irrelevant and we can 
consider only circularly-polarized face-on binaries. We 
fix the spin parameter to zero. 

Finally, as a reference we use the standard inspiral- 
only post-Newtonian approximation corresponding to the 
LAL approximant TaylorT4, which includes corrections 
to the phase of the waveform at 3.5PN order [TB] . 

The TaylorT4 waveforms used here are computed 
in the so-called "restricted" amplitude approximation, 
which assumes the waveform amplitude to be zeroth post- 
Newtonian order and only includes the leading second 
harmonic of the orbital phase. We only include the lead- 
ing second harmonic of the orbital phase in the EOBNR 
waveforms. We do not consider the effects of spin or ec- 
centricity in any of the waveform families, as we restrict 
to circular orbits and non-spinning black holes. 

The HG and Taylor T4 families are inspiral-only time- 
domain waveforms and are terminated when the gravita- 
tional waveforms reach the IS CO frequency. 



III. SNR FROM INSPIRAL, MERGER AND 
RINGDOWN 

In this section we consider the relative contributions 
of the different portions of the gravitational- wave coales- 
cence signal to the SNR as a function of the IMRAC's 
mass. 

We work in the frequency domain and define the 
Fourier transform of the gravitational- wave strain signal, 
h{f), as 



Hf) 



(3.1) 



where h{t) is the time-domain strain signal. We define 
the noise-weighted inner product as 



(a|6) = 43? 



^-^^ a{f)b*{f) 

/min ^nif) 



(3.2) 



where Sn{f) is the instrument noise power spectral den- 
sity (PSD), which we will take to be the Advanced LIGO 
high-power, zero-detuned noise PSD [20 . The limits of 
integration correspond to the bandwidth of the detector. 
The expectation value of the optimal matched filtering 
SNR, in the case when the signal and template wave- 
forms are identical, is given by ^I] 



/( 



f\Hf)\ 



d\nf 



1/2 



. (3.3) 



Writing the maximum SNR in the form above clearly 
separates it into contributions from the signal strain, 
f\h{f)\^ and the root-mean-squared (rms) noise spectral 
amplitude, \/ fSn^f)-, which is the strain signal associ- 
ated with the detector noise. 

One can gain insight into the relative contributions to 
the SNR from inspiral, merger and ringdown by com- 
paring the gravitational-wave strain to the noise rms 
value. In Fig. [T] we show the strain for selected overhead 
and face-on (i.e., optimally- located and oriented) IMRAC 
sources at a fiducial distance of IGpc as described by the 
EOBNR waveform family, and noise rms amplitude. The 
IS CO frequency of each signal is shown as a vertical line. 
The strain from merger and ringdown is thus the por- 
tion after the ISCO frequency. The contribution to the 
strain from merger and ringdown from binaries with com- 
ponent masses (mi, 7722) = [(200, 20) M©, (200, 2) M©], is 
greater than that of the noise rms amplitude (black curve 
with triangles). In general, systems with ISCO frequen- 
cies between 30 — 100 Hz, merge in the "bucket" of the 
noise curve, i.e., where the detector is most sensitive. 
For example, for the (mi, 7712) = (200, 20) system 
(red dotted curve in Fig. [T]), the merger and ringdown 
contribute the bulk of the SNR. Conversely, the inspiral 
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FIG. 1: Strain of optimally- located and oriented IMRAC sources at a fiducial distance of 1 Gpc as described by EOBNR 
waveforms, and aLIGO noise. The corresponding ISCO frequency of each signal is shown as a vertical dashed line. The strain 
from the merger and ringdown from each source contributes after the ISCO frequency. For the sources with component masses 
(7711,7712) = [(200,20) M0, (200,2) M0], the strain from merger and ringdown sits above the noise spectrum. The SNR from 
the full EOBNR waveform and from its inspiral-only portion are shown in Fig. [2] 



contribution to the SNR is strongly suppressed for such 
massive systems. 



In Fig. [2] we show the maximum SNR, Eq. (3.3), as a 
function of the binary's total mass produced by inspiral- 
only and full EOBNR waveforms at four different mass- 
ratios in the range 1/200 < q < 1/10. We construct 
inspiral-only EOBNR waveforms by Fourier transforming 
the full waveform into the frequency domain and trun- 
cating it at the ISCO frequency. We have considered 
the SNR for optimally-located and oriented sources at 
a fiducial distance of 1 Gpc. The lower-bound mass of 
the smaller body is set to 1712 = 1.4 M© which is the 
canonical neutron-star mass. The lowest total mass for 
the q = 1/50, 1/100 and 1/200 subplots is set by fixing 
the mass of the smaller body to 1712 = 1.4 Mq. For the 
q = 1/10 subplot in Fig. |2] the smallest total mass is set 
to M = 35 M0 as the inspiral phase accounts for the vast 
majority of the SNR below this value. The lower limit 



of integration of Eq. (3.3) is 10 Hz and the upper limit 
is 2048 Hz, which is the Nyquist frequency of discretely 
sampled EOBNR waveforms generated at a sampling rate 
of At = 4096s in the time-domain. We only consider 
systems with total masses such that the ISCO frequency 
is greater than 10 Hz (our low frequency cut-off). The 
highest total mass for each of the subplots in Fig [2] is 
set to M = 300 M0 which ensures the ISCO frequency 
is greater than 10 Hz. 

As anticipated from Fig. [l] there is a significant differ- 
ence in the SNR between inspiral-only and full EOBNR 
waveforms that can be seen at all four mass-ratios. We 
also note that for systems with mass-ratios of q = 1/10 
with total masses below around M = 35 Mq the inspi- 
ral phase is the dominant source of SNR. If we consider 
3% as a fiducial value of the difference between the full 
SNR and the one associated to the inspiral-only wave- 
form - which leads to a loss in detection rates of 10% - 



5 



this happens ai M ^ 35 Mq for q = 1/10. For binaries 
with q = 1/50, 1/100 and 1/200, the minimum differ- 
ence in SNR between inspiral-only and full waveforms is 
~ 6%, 15% and 40% respectively for the mass ranges 
considered in Fig|2] 

In summary, we have shown that inspiral-only tem- 
plates will miss a significant portion of the total SNR 
of IMRAC signals over the bulk of the detectable mass- 
space. Future searches will therefore require templates 
that can match the full inspiral-merger-ringdown. How- 
ever, there is a small region of the parameter space for 
which inspiral-only templates may suffice for searches, 
without inducing drastic losses in detection rates. In the 
following section we quantify the effectiveness of inspiral- 
only templates for searching for full coalescence signals 
in aLIGO. 



IV. EFFECTIVENESS OF INSPIRAL-ONLY 
TEMPLATES FOR IMRAC SEARCHES 



We have shown in Sec. Ill that the SNR from merger 
and ringdown will provide a significant contribution to 
the total SNR over a broad portion of the IMRAC mass- 
space, c.f. Fig. |2] There is however a small portion of 
the parameter space where the SNR is dominated by the 
inspiral phase. This can be seen in Fig. |2] for binaries 
with q = 1/10 binaries in Fig. |2] with total masses M < 
35 Mq . Thus it is important to quantify the effect of us- 
ing inspiral-only templates to search for IMRAC signals 
which contain inspiral, merger and ringdown phases. The 
use of template waveforms that are not exact representa- 
tions of the signals they filter degrades the SNR, as the 
optimal SNR can be recovered only when the template 
waveform corresponds exactly to /i, see . Eq. (3.3). In 



practice however, we do not have access to an exact rep- 
resentation for h. Using a non-exact template waveform 
T to filter h caps the maximum recoverable SNR to 



max - 



iH\)\m) 

(T(^)|T(^))i/2 

max 



(4.1) 



where A and represent the parameter vector of the sig- 
nal and template, respectively. We define e as the effec- 
tiveness of a template waveform family T at recovering 
the maximum SNR from a gravitational- wave signal h; 
by definition < e < 1. This quantity is also referred 
to as the "fitting factor" in the literature [22]. It is con- 
venient to define waveforms normalized to unit norm as 
= a(/)/(a|a)i/2 so that (h|h) = (t|t) = 1 and the 
effectiveness can be written succinctly as [22^ 



max(h(A)|t((9)). 



(4.2) 



Using normalized waveforms also has the advantage of 
eliminating the dependence of the waveforms on the 



source orientation and distance, which enter as an overall 
scaling. 



Calculating the effectiveness, Eq. (4.2), requires maxi- 



mizing over the component masses (mi, 777.2) and the time 
and phase at coalescence. We can efficiently maximize 
over the time and phase by Fourier transforming the in- 
tegrand of the noise- weighted inner-product ^23J, 



Z{tc) 



frr 



frr 



^ Snif) 



(4.3) 



which yields a complex time series whose elements cor- 
respond to the inner-product of a and b as one of the 
signals is time-shifted with respect to the other. We can 
efficiently find the time at coalescence, tc, by finding the 
time at which the norm of this time series is a maximum. 
The phase at coalescence <pc is then automatically given 
by finding the argument of the time-series at its peak 
amplitude. We thus modify the inner product {a\b): 



{a\b) ^ {a\by 



niax \z{tc 



(4.4) 



which we will adopt as the definition of the inner-product 
for the remainder of this paper. 

To compute the effectiveness of an inspiral-only IM- 
RAC search we evaluate Eq. (4.2) for signals covering 



the IMRAC mass space. We take as our signal waveform, 
/i, the full inspiral-merger-ringdown EOBNR waveform. 
We take the template, T, to be an inspiral-only EOBNR 
waveform, formed by truncating the full EOBNR wave- 
form at the ISCO frequency in the frequency domain. 
With such signals and templates the effectiveness pro- 
vides a measure of the maximum SNR which could be 
achieved through using an inspiral-only template to fil- 
ter full coalescence-signals. To get a br oad coverage of 
the IMRAC mass space we compute Eq. ( |4.2| ) for signals 
whose source masses cover the ranges 1.4 Mq < m2 < 
18.5 M© and 24 M© < mi < 200 M©, with mass ra- 
tios spanning the range q := m2/mi G [1/140,1/10]. 
For each signal we evaluate the effectiveness, Eq. ( |4.2[ ), 
where the template T describes the inspiral-only por- 
tion of an EOBNR waveform. We maximize over time 
and phase by maximizing the inner product of the signal 
with an inspiral-only EOBNR template, Eq. (Ol). The 



maximization over the masses is performed by finding 
the largest inner product between the signal and a bank 
of template waveforms constructed so that the minimal 
inner product between (normalized) neighbouring tem- 
plates is 99%. The template bank is characterised by in- 
trinsic parameters (mi, m2) and spans an extended mass 
range mi ^2 x (1 i 0.1) where mi ^2 are the masses as- 
sociated with each signal waveform. The results of the 
effectiveness of an inspiral-only IMRAC search are shown 
in Fig. [3| 

Inspiral-only templates are ~ 98% effective at filter- 
ing full coalescence signals for total masses M ^ 30 Mq. 
Such systems have an ISCO frequency 150 Hz ^ /isco 
which is well within the peak sensitivity of the noise 
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FIG. 2: SNR of optimally-located and oriented IMRAC sources at a fiducial distance of 1 Gpc vs total mass for four different 
mass ratios; q — 1/10, 1/50, 1/100, 1/200. The solid line is the SNR from full EOBNR waveforms and the dashed line from 
inspiral-only EOBNR waveforms truncated at the ISCO frequency in the frequency domain. The lower-bound mass of the 
smaller body is set to m2 — 1.4 M© which is the canonical neutron-star mass. The lowest total mass for the q — 1/50, 1/100 
and 1/200 subplots in Fig. [5] is set by fixing the mass of the smaller body to m2 = 1.4 M©. For the q — 1/10 subplot the 
smallest total mass is set to M = 35 M© as the inspiral accounts for the vast majority of the SNR below this value. We only 
consider systems with total masses such that the ISCO frequency is greater than 10 Hz (our low frequency cut-off). The highest 
total mass for each of the subplots in Fig |2] is set to M = 300 M© which ensures that the ISCO frequency is greater than 
10 Hz. We find that there is a non-negligible contribution to the SNR from merger and ringdown in IMRAC signals above a 
total mass of around M = 35 M© . The difference in SNR between inspiral-only and full waveforms is at the 3% level at around 
M = 35 Mo at = 1/10. For binaries with q = 1/50, 1/100 and 1/200, the minimum loss in SNR are at the 6%, 15% and 40% 
levels, respectively, in our mass range of interest. For IMRACs of astrophysical interest, more extreme mass ratios correspond 
to greater total mass, which can place the merger and ringdown at a frequency where the detector has the greatest sensitivity. 




curve. However for the bulk of the mass space the ef- 
fectiveness is below 75%. This is unsurprising given the 
SNR curves in Fig. |2] which clearly show the importance 
of the contribution of merger and ringdown to the SNR. 

The loss in SNR incurred through using inspiral-only 
templates directly affects detection rates. Because the 
SNR scales inversely with the distance, the observable 
volume will scale with the cube of the effectiveness. As- 
suming that GW sources are isotropically distributed in 
the sky, the fractional loss in detection rates will be 
1 — e^. The percentage loss in detection rates through 
using inspiral-only EOBNR templates to recover the full 
coalescence signal is also shown in Fig. [3| Over a broad 
portion of the mass-space inspiral-only templates incur 
losses in detection rates between 60 — 85%. As the total 
mass of the binary approaches 440 the Schwarzschild 



ISCO frequency, Eq. ( |1.1| ), approaches 10 Hz which is 
near the low-frequency cut-off of the detectors. Hence 
the relative contribution of the inspiral phase to the co- 
alescence signal of heavier systems diminishes until the 
only contribution is from the merger and ringdown. This 
is a striking indication of the need of merger and ring- 
down in IMRAC template waveforms. This suggests the 
importance of full numerical simulations in this regime in 
order to construct a reliable waveform family including 
inspiral, merger, and ringdown phases. 

We identify three regions in the mi-m2 plane in which 
various searches could be constructed. The regions are 
defined by contours of constant effectiv eness wh ich are 
approximately given by C = {mi / Mq) ^/m2JM^ ^ which 
are found purely empirically, with 1.4 Mq < m2 < 
18.5 M0 and mass-ratios q G [1/140, 1/10]. The effective- 
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FIG. 3: Effectiveness of inspiral-only EOBNR templates to filter full inspiral, merger and ringdown EOBNR signals as a 
function of the source component masses and corresponding loss in detection rates. The diagonal corresponds to a mass-ratio 
q = 1/10. Inspiral-only EOBNR templates are constructed by truncating the full waveform at the ISCO frequency in the 
frequency-domain. For the bulk of the parameter space inspiral-only templates are ^ 75% effective at filtering inspiral, merger 
and ringdown signals. Inspiral-only templates are ^ 97 — 98% effective for total masses Mq ^ 30 Mq. Inspiral-only templates 
within the 90%-effectiveness contour should be sufficient for IMRAC searches without incurring greater than 30% losses in 
detection rates. 



ness is related to C by e 1/100 x (1.6 C- 7.3 x 10"^ C^). 

Between the 97%- and 90%-effectiveness contours, the 
losses in detection rates are between 10% ^ L ^ 27% 
and so an inspiral-only search could be sufficient without 
incurring drastic losses in detections. The region bound 
from below in effectiveness by the 90%-effectiveness con- 
tour is defined by C < 100, with the effectiveness in- 
creasing with decreasing C. Between the 90%- and 80%- 
effectiveness contours the losses in detection rates are 
around 27% ^ L ^ 50%. Thus, within this region 
searches will be limited by the lack of merger and ring- 
down in template waveforms, though an inspiral-only 
search would be feasible in principle. This contour is 
defined by 100 < C < 150. Below the 80%-effectiveness 
contour, inspiral-only searches will incur losses in detec- 
tion 50% < L and so merger and ringdown will be crucial 
for searches. The region bound from above in effective- 
ness by the 80%-effectiveness contour which is defined by 
150 ^ C, with effectiveness decreasing with increasing C. 
The results are summarized in Table [H 



V. COMPARISON OF INSPIRAL-ONLY 
WAVEFORMS 

We have shown that merger and ringdown are crucial 
for effective searches over a large portion of the IMRAC 



mass space, though there is a small region in which an 
inspiral-only search could be constructed without incur- 
ring losses in detection rates greater than around 27%. 
For this region, it is therefore important to study whether 
currently available waveforms are sufficiently accurate. 
The inspiral phase can be computed using perturbative 
expansions and thus it is interesting to quantify the con- 
sistency of different expansions. 

To assess the effectiveness of the EOBNR inspiral, 
we employ a waveform family designed to approximate 
intermediate mass-ratio inspirals which we refer to as 
"Huerta-Gair" (HG) waveforms [15 . HG waveforms de- 
scribe only the inspiral portion of the coalescence signal. 
This waveform family has no corresponding LAL approx- 
imant. 

We repeat the study done in the previous Section us- 
ing now the HG waveform family as the signal h and 
inspiral-only EOBNR as the template T. The results are 
reported over the whole parameter space in Fig. [4] For 
completeness, in Table [ll] we also show the values of the 
effectiveness, Eq. |4.2[ for selected mass combinations of 
EOBNR inspiral-only templates for filtering full EOBNR 
and HG signals respectively. 

For the high-mass part of the mass-space the effec- 
tiveness of the EOBNR inspirals with respect to the HG 
waveforms is close to 100%. This is perhaps unsurprising 
because very high mass systems will have short inspi- 
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Effectiveness of 
inspiral-only search, 


Loss in detection 
rates, L{%) 


Contours of constant effective- 
ness in mi -7712 plane {C — 
[mi/ MQ)^m2/ Mq) witnm mass 
range of interest 


Implication for searches 


90% < e < 97% 


10% <L< 27% 


C < 100 


Inspiral-only search sufficient but 
with non-negligible loss in detection 
rates. 


80% < e < 90% 


27% < L < 50% 


100 < C < 150 


Inspiral-only search possible but 
limited by lack of merger and ring- 
down. Could potentially lose half 
of signals with inspiral-only tem- 
plates. 


e < 80% 


50% < L 


150 <C 


Merger and ringdown crucial for 
searches. Will miss over half of sig- 
nals with inspiral-only templates. 



TABLE I: Effectiveness of inspiral-only searches, the corresponding loss in detection rates and the region in the mi-m2 plane 
bounded by constant-effectiveness contours. For a given region in the mi-m2 plane bounded by const ant- effectiveness contours 
we summarize the implications for IMRAC searches. 
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FIG. 4: Effectiveness of inspiral-only EOBNR templates at filtering HG signals (left) and TaylorT4 signals (right) as a function 
of the source component masses encoded in the signal. 



rals and possible differences in the waveforms will not 
produce a significant degradation of SNR when matched 
over a small number of wave cycles. However, for lighter 
systems the effectiveness can be as low as 90%, which 
occurs in the region of mass space in which inspiral-only 
searches would be most feasible (see Table 

For reference we also compare inspiral-only EOBNR 
templates to TaylorT4 signal waveforms (which are 
inspiral-only). We construct signal waveforms on the 
same grid in mi — m2 as for HG waveforms and use the 
same template bank of inspiral-only EOBNR waveforms. 
The results are summarized in the right panel of Fig. |4j 
and in Table [ill for selected masses. We find that the 
EOBNR inspiral has good filtering efficiency for the Tay- 
lorT4 waveform family. However, EOBNR is clearly a 
better match to the HG waveform family over a larger 
range of masses and mass ratios than to TaylorT4. This 



can be seen more clearly by comparing the subplots in 
Fig. [4] This is unsurprising given that the PN expansion 
is unreliable at high velocities and highly asymmetrical 
mass-ratios. For orbital velocities v/c = (Mtt/)^/^ ^ 0.2 
the PN energy flux deviates significantly from numerical 
results, see [e.g.,|24l[2S)- A binary at its ISCO frequency 
has v/c ^ 0.4, which is well beyond the region of validity. 



VI. DISCUSSION AND CONCLUSION 

We have shown that over the bulk of the IMRAC mass 
space, merger and ringdown contribute significantly to 
the gravitational-wave coalescence signal. This happens 
despite the suppression of the power in the merger and 
ringdown in signals from binaries with very asymmetric 
mass ratios. The importance of merger and ringdown is 
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mi 


m2 




Signal waveforms 




(Mo) 


(Mo) 


full EOBNR Huerta-Gair TaylorT4 


50 


5 


0.90 


0.95 


0.96 


100 


5 


0.76 


0.97 


0.97 


200 


5 


0.53 


0.99 


0.99 


50 


1.4 


0.96 


0.96 


0.90 


100 


1.4 


0.86 


0.98 


0.94 


200 


1.4 


0.67 


0.98 


0.98 


5 


5 


0.99 


0.89 


0.99 


20 


20 


0.99 


0.92 


0.99 


100 


100 


0.52 


0.51 


0.50 



TABLE II: Summary of effectiveness of inspiral-only EOBNR 
template waveforms in recovering signals modelled using dif- 
ferent waveform families - full EOBNR, Huerta-Gair and Tay- 
lorT4 - for selected component masses. Merger and ringdown 
become more prominent in the coalescence signal as the total- 
mass of the system is increased. The EOBNR inspiral is typ- 
ically better at matching HG signals in the IMRAC regime 
than TaylorT4 signals. Results for equal mass-ratio systems 
are shown for reference below the horizontal line. 



and a region in which merger and ringdown are critical 
to prevent losses in detection rates over 50%. The search 
regions are summarized in Table [Tj 

We have further shown that in the region of the IM- 
RAC mass space in which inspiral-only searches are fea- 
sible, approximants adapted to asymmetric mass-ratio 
binaries are important, as here the binary is liable to 
have highly relativistic velocities v/c ^ 0.2. We consid- 
ered a waveform family designed to describe intermediate 
mass-ratio binaries which we referred to as the "Huerta- 
Gair" (HG) waveform family. By computing the effec- 
tiveness of inspiral-only EOBNR waveforms to filter sig- 
nals described by the HG waveform family, we showed 
that losses in recovered SNR could be as great as 10%. 
In Table [H] we summarize the effectiveness of the signal- 
template combinations used in the paper. 

We believe that template waveforms for IMRAC 
searches will benefit from calibration to several numerical 
simulations. We note that there already exists one very 
short numerical waveform oi a q = 1/100 binary which 
we have not used in our study, and which EOBNR is not 
currently calibrated to [llj . 



due to the greater sensitivity to these waveform portions 
for high-mass signals, for which most of the inspiral may 
fall at frequencies below the detector's sensitive band. 

However there is a relatively large patch in mass space 
in which the inspiral-only waveforms are more than 90% 
effective. We identified three regions in which differ- 
ent searches could be considered appropriate based on 
thresholds of acceptable losses in detection rates. The 
mass space splits into a region in which inspiral-only 
searches could be feasible, incurring losses in detection 
rates of up to ~ 27%; a region in which searches would 
be limited by lack of merger and ringdown in template 
waveforms, incurring losses in detection rates up to 50%; 
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